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Abstract 

We apply the formalism of quantum cosmology to models containing a phantom field. Three 
models are discussed explicitly: a toy model, a model with an exponential phantom potential, and 
a model with phantom field accompanied by a negative cosmological constant. In all these cases we 
calculate the classical trajectories in configuration space and give solutions to the Wheeler-DeWitt 
equation in quantum cosmology. In the cases of the toy model and the model with exponential 
potential we are able to solve the Wheeler-DeWitt equation exactly. For comparison, we also give 
the corresponding solutions for an ordinary scalar field. We discuss in particular the behaviour of 
wave packets in minisuperspace. For the phantom field these packets disperse in the region that 
corresponds to the Big Rip singularity. This thus constitutes a genuine quantum region at large 
scales, described by a regular solution of the Wheeler-DeWitt equation. For the ordinary scalar 
field, the Big-Bang singularity is avoided. Some remarks on the arrow of time in phantom models 
as well as on the relation of phantom models to loop quantum cosmology are given. 
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I. INTRODUCTION 



It is a striking fact that our universe is currently accelerating A major open problem 
is to provide a fundamental type of matter which may be responsible for this, since any of 
the forms of matter we know from our experience cannot explain this phenomenon. This 
matter is not visible, but provides a dominant fraction of the energy density in the universe 
and was therefore given the name 'dark energy' (for a current review see for example 
The best known, and perhaps the simplest candidate for such a matter is the cosmological 
constant, but theoretical physics provides more options. One of them is an evolving scalar 
field with appropriate kinetic and potential energies. In general, it may mimic various types 
of matter during different periods of the cosmological evolution. 

Dark energy is characterized by negative pressure which causes the repulsion of matter 
in the universe and, as a consequence, its acceleration. In terms of the standard energy 
conditions in general relativity |3|, dark energy must violate the strong energy condition 
p + 3j9 > 0, p > 0. Assuming a barotropic equation of state of the matter in the universe, 
p = wp {w = constant), where p and p are the pressure and the density of dark energy, 
respectively, it requires that w < — |. 

However, according to more recent observations 0, |^ , dark energy is even more biased 
towards larger negative values of the barotropic index w < — 1. This means that it would 
have to violate the null energy condition g + p > and, consequently, all the remaining 
energy conditions such as: the weak energy condition g > 0, g + p > 0, and the dominant 
energy condition g > 0, —g < p < g. Dark energy of this type was dubbed phantom 0|. 
A phantom may be represented by an evolving scalar field which possesses negative kinetic 
energy (often called 'ghost field'). Although phantom fields lead to various problems j^, 
they are observationallv supported as a possible source of dark energy and deserve thorough 
investigation (but see [2| for an alternative view). Moreover, there may exist phantom fields 
without pathologic behaviour in the ultraviolet regime (lo| . 

. ^hantom models of the universe admit a new type of singularity called a Big- Rip singularity 
0, llH- At the Big Rip, energy density and pressure diverge as the scale factor a{t) 
goes to infinity at a finite time. This is different from an ordinary Big Crunch singularity, 
which leads to a blow-up of the energy density and pressure as the scale factor approaches 
zero at a finite time. Another possible singularity is the Big Brake where the expansion 
rate is zero and the acceleration rate approaches minus infinity Besides, more exotic 
types of singularities may appear such as the sudden future singularity [l^, the generalized 
sudden future singularity [l^ where there is a blow-up of the higher-order derivatives of the 
scale factor with smooth evolution of the scale factor and the energy density, the type III 
singularity and the type IV singularity where the evolution of the scale factor is 
smooth. These singularities have weaker properties than a Big Rip [l7| . 

In Ref. classical phantom cosmologies were studied and a large variety of possible 
cosmological scenarios were found. Also, the dualitv relation between standard matter and 
phantom matter models was revealed (see also_|l9|,|20|) which has an analogon in the duality 
symmetry present in superstring cosmology [2l|. 

It is worth mentioning that, once the supernova data are analyzed in a prior-free manner. 
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an evolving equation of state with a time-dependent barotropic index w = w{t) for the dark 
energy is favoured (cf. [22, |23|, |2J]). Such models were also studied in the quantum context 
in [25| where a canonical momentum was attached to a time-dependent barotropic index. 
Some authors have also studied the thermodynamical properties of phantom models ji^ li^l . 

In all these investigations, an evolving universe was described by classical cosmology. 
Quantum effects were only studied in certain phases of the evolution, for example, close to 
a singularity, cf. ji^, without applying quantum theory to the universe as a whole. It is the 
purpose of this paper to accomplish the latter goal, that is, to discuss quantum cosmology 
with phantom fields. The interest in this is due to the fact that for both experimental 
and theoretical reasons it seems that quantum theory is universally valid [2^. Therefore, 
the universe as a whole has to be described by quantum theory. If phantom fields play a 
dominant role, it has to be investigated whether this causes deviations from the standard 
formalism of quantum cosmology and whether there are interesting physical consequences. 

Quantum cosmology must be based on a theory of quantum gravity [s^ . Candidates for 
such a theory include string theory, loop quantum gravity and quantum geometrodynamics. 
Our present analysis will, like most investigations of quantum cosmology, be based on the 
Wheeler-DeWitt equation of quantum geometrodynamics. Independently of the correct 
theory of quantum gravity, this framework should yield an adequate description at least on 
the energy scales below the Planck scale (if not on all scales) . If one approaches the Planck 
scale, modifications such as loop quantum cosmology [sH might become necessary. The 
investigations in our paper are independent of such modifications and will be discussed in a 
future paper. 

A central feature of the Wheeler-DeWitt equation is its local hyperbolic signature 
In regions of configuration space near closed Friedmann cosmologies, it is globally hyperbolic, 
that is, there is only one minus sign in the kinetic term [2^. The negative part of the kinetic 
term is related to the scale factor of the Friedmann model, which in a certain sense thus plays 
itself the role of a phantom field. The presence of an indefinite kinetic term is intimately 
connected with the attractive nature of gravity [j^ ]. 

Besides its hyperbolic character, the most important feature of the Wheeler-DeWitt 
equation is its independence of an external time parameter 3^ S^l- This holds, in fact, for 
every system that is reparametrization invariant at the classical level. Consistent discussions 
of quantum cosmology must thus be based on the intrinsic structure of this equation and 
avoid the use of an intuitive but wrong picture of an external Newtonian time. For this 
purpose it is necessary to study the classical trajectories in a configuration space where the 
classical time parameter t is eliminated. 

The structure of the Wheeler-DeWitt equation is important for the imposition of bound- 
ary conditions in quantum cosmology. In the hyperbolic case one has a wave equation whose 
form suggests imposing boundary conditions at constant values of the scale factor, a. This 
is of importance, for example, if one attempts to construct wave packets that follow the clas- 
sical trajectories in configuration space like standing tubes jsi, 3^, 3^. It is also crucial for 
an understanding of what pre- and post big bang phases mean in quantum string cosmology 
Isii liil. The origin of the arrow of time can in principle be traced back to the structure of 



this wave equation 32, 4^ 
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The presence of a phantom field changes the structure of the Wheeler-DeWitt equation: 
If only the phantom is present besides the scale factor ('phantom dominance'), its structure 
becomes elliptic, while in the general case it becomes of a mixed (ultrahyperbolic) nature. 
This has implications for the imposition of boundary conditions. A change of signature in 
the Wheeler-DeWitt equation has hitherto been noticed in the presence of non-minimally 
coupled fields i^]- In our paper we shall present the formalism of quantum phantom cos- 



mology and some of its main physical consequences. 

Our paper is organized as follows. In Sec. |n] we shall study and solve the classical 
equations of motion for the phantom field in a Friedmann universe. After the presentation 
of the necessary equations, we give the solutions for the classical trajectories in configuration 
space for three models: A toy model with vanishing phantom potential (Sec. IIIB|1 . a model 
with exponential phantom potential (Sec. Ill C|) . and a model with cosh-potential and a 
negative cosmological constant (Sec. IIID|) . For comparison, in all these cases, we give the 
results for a non-phantom scalar field. Sec. IIIII contains in the same order the discussion of 
the quantum theory for these models, both for a phantom field and a corresponding ordinary 
scalar field. We are able to solve the Wheeler-DeWitt equation exactly for the toy model 
and the model with exponential potential. In particular, we discuss wave packet solutions 
and find that quantum effects dominate in the region of the classical Big-Rip singularity. 
Therefore, quantum effects occur at large scales. Since the solutions of the Wheeler-DeWitt 
equation are regular there, the Big-Rip singularity has vanished in the quantum theory. 
Furthermore, in the realistic scalar field models, the wave function vanishes at the Big 
Bang. Thus, this singularity is likewise excluded in the quantum theory. In Sec. HXlwe give 
a summary of the results and the outlook of the problem of the arrow of time and possible 
modifications of the obtained picture due to loop quantum cosmology. 



II. CLASSICAL PHANTOM COSMOLOGIES IN CONFIGURATION SPACE 
A. Classical equations 

We consider a Friedmann universe with scale factor a{t) and a homogeneous scalar field 
(t){t). We assume here that the 0-field dominates over other matter degrees of freedom, so 
that it is the only degree of freedom besides the scale factor. The action reads 

Here, = 8ttG, N is the lapse function, A is the cosmological constant, V{(j)) is a potential 
of the field 0, /C = 0, ±1 is the curvature index, and we have set c = 1. The parameter i 
distinguishes between a phantom field (where i = —1) and an ordinary scalar field (where 
i = +l). 

We set = 1, so the time parameter is the standard Friedmann cosmic time. The action 
then becomes 

^ = A f dt {-aa^ + lCa- ^a^) + \ ! dt {a^i<i)^ - 2aV((/>)) . (2) 
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The canonical momenta are given by 

Qad 



TT^ = ia'^<j) . (3) 



The canonical Hamiltonian Ti, which is constrained to vanish, reads 

K.^ 2 ^'^l 3/Ca ^ 

n = -—7Tl + -^ + a^— + a^V ^ = 0. 4 

12a 2a^ h? 

Expressed in terms of the 'velocities', see (jH)), this constraint becomes identical to the Fried- 
mann equation. 

The term in parentheses is the energy density of the scalar field. 

We recognize that for the standard scalar field (£ = 1), no classically forbidden regions in 
configuration space exist due to the indefiniteness of the total kinetic term. This is different 
from the phantom case {^ = —1), where only the region 

is classically allowed (this restriction is due to the negative definiteness of the total kinetic 
term). 

The field (p obeys the second-order equation of motion 

+ 3-0 + £V^'(0) = , (8) 
d 

which is equivalent to the conservation equation p + 3H{p + p) = 0, provided the standard 
perfect-fluid energy-momentum tensor is introduced. This equation is trivially fulfilled by 
the cosmological constant A (p = constant, p = — p). In (jHJ we recognize a formal reversal 
of the potential in the phantom case compared to an ordinary scalar field case, since the 
sign in front of the V^'-term changes. With the help of ((HI), the second-order equation for a 
can be put into the form 

^t4(^^^-VW)^0. (9) 

Again, assuming the perfect-fiuid energy-momentum tensor, the scalar field exerts the pres- 
sure 

p^i^- V{<jy) . (10) 

Note that the case of a cosmological constant is included by having the additional equation 
of state 

A , , 

Pa = -Pa = ^ • (11) 
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Assuming a constant barotropic index w for the scalar field, we can use (jHl) and (fTUj) to find 
the relation between the scalar field and its potential 43[, 



Vim 



w 



21 + w 



(12) 



This is analogous to the virial theorem in which the kinetic energy is proportional to the 
potential energy of the field. However, as has already been mentioned above, it may be more 
physical to assume a time- dependent barotropic index 22, |23|, |24 



B. Classical phantom trajectory for vanishing phantom potential and vanishing 
cosmological constant 

In this section we shall consider a simple model with field potential V{(j)) = and cosmo- 
logical constant A = 0. This leads to an equation of state for stiff matter, p = p, w = 1, in 
contrast to the current observational status jl[ 1^ ■ However, such an evolution may perhaps 
be valid in ekpyrotic/cyclic scenarios where this matter dominates the collapsing phase of 
the cosmological evolution Q|. Moreover, in such a case the energy density p < 0, and 
thus this model does not seem to represent dark energy which is usually assumed to have 
positive energy density. However, it captures interesting 'phantom features', since it violates 
all energy conditions, and it has the merit that it is easily manageable. More realistic models 
will be discussed below. 

In order to get classical solutions in the phantom case [i = —1) we have to choose 
/C = — 1 in (jZj). Since we are interested in constructing wave packets from the Wheeler- 
DeWitt equation, we want to find a classical trajectory in configuration space, where the 
classical time t is eliminated. This is motivated by the fact that no such parameter is present 
in the Wheeler-DeWitt equation. 

Since is a cyclic variable, tt^ is constant, so from one has 0^ = C^/a^ with a constant 
Cp. From (0) one then has (we choose Cp > 0) 



dcj) , C, 



p 



(13) 



a" - 



which can easily be integrated to yield 



0(a) = i-W^arccos-^^ft; • (1^) 



V6i 

For convenience, we choose = 6. Then the solution reads 

1 C 

(l)(a) = ±-arccos^ . (15) 

2 a'' 

The classical trajectory (fT^ has a minimum value of the scale factor, amin = a/C^? and 
reaches infinite values of a at finite values of = ±7r/4. In this sense it resembles a Big- Rip 
solution. However, with respect to t the scale factor reaches infinity only at t = ±oo and, 
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FIG. 1: The classical trajectory in configuration space for the toy model with vanishing scalar 
field potential and vanishing cosmological constant. The diagram on the left-hand side shows the 
trajectory for the phantom field model. On the right-hand side the trajectory for the i = 1 scalar 
field model is plotted. 



moreover, p oc which is the density scaling appropriate to a stiff fluid. Nonetheless, in 
configuration space the trajectory has some features of a Big-Rip, and this is why this toy 
model is of interest. 

For an ordinary scalar field {£ = 1) and for /C = — 1, one gets instead of (fT3jl . 

1 C 

0(a) = ±-arcsinh— ^ . (16) 

2 ci 

There is no turning point; Eq. just describes two branches for which a — oo if — > 0, 
and a —>■ if (p ±oo. The two solutions (fTH|) and pfij) are depicted in Figure 1. For /C = 1 
one obtains the solution with a turning point (arccosh instead of arcsinh) that was discussed 
in 



C. Classical trajectories for exponential scalar field potential and vanishing cos- 
mological constant 

A model with phantom equation of state and a true Big-Rip singularity for the phantom 
model appears if the potential in ^ is chosen to be exponential 



45, 4 



V (0) = Voe-^'^^ , (17) 

and A = 0. Interest in this type of scalar field potentials in cosmology arose when it 
)ecame clear that the classical model has an attractor solution with scalar field domination 
This alleviates the fine-tuning problem of the initial energy of the scalar field 



48l |. Such an attractor exists not only in the case of a conventional scalar field, but also 
for the phantom field 0, 4^|- Exponential potentials for scalar fields arise in the context 



of Kaluza-Klein theories [sCll EH , higher-derivative gravity in {D + 4) dimensio ns [5^ 0| , 
higher-order gravity , supergravity and superstring theories Eo] , see also jijfor an 



overview. 
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In the following, we shall consider the case of a flat universe, = 0. From equation 
([7j) one sees immediately that for this choice of parameters neither the ordinary scalar field 
nor the phantom field model possesses classically forbidden regions. The classical equations 
of motion (jHJ and Q can be transformed into a dynamical system with the Friedmann 
equation © constraint 0,0,143,13. For i = —1 and arbitrary values of A, as well 
as for i = +1 and A < y/E, this system has an attr actor solution given by [l9l |: 



1 + 



(t - to) 



(18) 



a 
ao 



1 + 



(19) 



Introducing a = ln(a) for later convenience, one obtains the following simple trajectory in 
configuration space, 

(j){a)=£-a. (20) 



For this attractor solution, the 'kinetic energy' defined from (0) — writing this equation 
in the form £-Ekin + -Spot = 1 — is given by (using = 6 in the second step) 



kin 



~6 



da 



A^ 
6 



and thus constant. Therefore, also the 'potential energy' of the scalar field is constant, 

K^V 



pot 



1 



£A^ 
6 



The equation of state parameter of the field is 



w 



X2 



(21) 



Thus for 
field with 



1, indeed describes a phantom field with w < —1, whereas the scalar 
1 covers the range w > —1. Accordingly, the energy density scales as 

p = Po y-^j ■ As expected, this yields a Big- Rip singularity for i = —1, since in the limit 

t ^ ti = to — 2C,/ (X^Hq) the energy density and the scale factor diverge. For t — oo, a and 
p vanish. This is in contrast to the i = 1 model: In the limit t ti, a goes to zero and p 
diverges, yielding a Big Bang, whilst for t ^ oo, a diverges and p goes to zero. 



D. Classical trajectories for scalar field fluid and negative cosmological constant 

It is easy to obtain a simple set of classical solutions for cosmological models with a 
negative cosmological constant In contrast to a positive cosmological constant which 
supports cosmological repulsion, the negative cosmological constant is a source of attraction 
and can overcome the influence of repulsion from dark energy with negative pressure such 
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as cosmic strings, domain walls, and phantom, see for example [18[- This allows models 
with a negative cosmological constant and other fluids to evolve symmetrically between two 
singularities with an extremum in between. In particular, it is possible to have an evolution 
beetween the two Big Rips which appear at finite cosmic time, as will be shown below. 

We assume a fiat universe, /C = 0, with a negative cosmological constant A < and a 
cosmological fluid with barotropic equation of state p = wp; the latter will be mimicked by 
a scalar field (p of either standard or phantom type. In this case, the energy conservation 
equation gives 

p = Ca-''^''"+^'> , (22) 

and the energy conservation equation (|TT| for the cosmological constant remains valid. This 
can be used to solve the system (0) and Q in terms of the scale factor as 



a{t) 

where 



(23) 



Q an 

D = -{l+w), — = -A>0. (24) 
Using ()24|). we can rewrite ()12|) in the form 

Vim = ' (25) 

which allows to write (jH} as (note that as in Sec. Ill Bl we have assumed = 6), 

With all these assumptions we are able to calculate the evolution of the scalar field as 

0(t) = ±^i— \/ V In I tan ( ^{-A)n] \ . (27) 



Let us note that i = +1,D > for an ordinary scalar field, while i = -~1,D < for 
the phantom. Then, the above expressions make sense since D/i = \D\ > 0. For D > 
(negative cosmological term plus w > —1 fluid), the evolution of the universe based on 
begins with a Big Bang at t = 0, reaches a maximum flmax = A^^^, and terminates with 
a Big Crunch at t = vr. For D = —\D\ < (the phantom case), the evolution starts with 
a Big Rip at t = 0, reaches a minimum a^i^ = A"^/'^', and terminates with a Big Rip at 
t = IT. The latter case is of special interest because it allows a symmetric evolution of the 
scale factor in the presence of a phantom field. This model may also be of interest to study 
the cosmological arrow of time, see Sec. IV. 

A similar type of symmetric evolution appears in configuration space. Using ()23j) and 
()27|1 to eliminate the classical time coordinate, we obtain the trajectory 




a" 



^("'^^^^/^'■'IwF^l' '''' 
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FIG. 2: The classical trajectories in configuration space for the models with cosh-potential and 
negative cosmological constant. On the left-hand side, the trajectory for the phantom field model 
is shown. The classical trajectory for the scalar field model £ = 1 is shown on the right-hand side. 
The similarity to the classical trajectories in the toy model in Sec. IIIBl is obvious. 



From this we can see that there are two branches. For £ = —1 each of them extends to 
infinity, that is, (p — > ±oo, for a ^ oo and reaches a minimum 0(a) = 0, for amm = A~^^^^^. 
For i = 1 one recognizes the presence of the maximum Omax- The trajectories in configuration 
space are depicted in Figure 2. From ^I7\i and (pUj) one can reconstruct the potential of the 
scalar field, 

\/(0) = Hcosh2(|) ^ (29) 

where 

A / 1 

Note that for £ = 1 the potential is positive only for D < 3 (i.e., w <1). This restriction is 
similar to the restriction A < a/6 in Sec. Ill CI 




III. QUANTUM COSMOLOGY FOR PHANTOM AND ORDINARY FIELD 
A. Wheeler DeWitt equation and phantom duality 

Quantization of the Hamiltonian constraint ^ leads to the Wheeler-DeWitt equation. 
Choosing the Laplace-Beltrami factor ordering and again the convention = 6, it reads 

Let us note that under the phantom duality [l^ 

a ^ - , (31) 

a 

^ -20 , (32) 



10 



for /C = the Wheeler-DeWitt equation for a, 0, 

transforms into the Wheeler-DeWitt equation for a, 0, that is, 

The transformation for is thus just a Wick rotation. 

On the other hand, Eq. (jH^ can conveniently be rewritten in terms of the scale factor 
a = In (a) as 

(fl^ - 41^ + (^^) 

It is this form of the Wheeler-DeWitt equation with which we shall work in the following. 



B. Quantum phantom cosmology — no phantom potential 

For vanishing potential = 0, A = 0, and /C = — 1 the solution to the phantom {i = — 1) 
Wheeler-DeWitt equation (|Hn|) is found by a separation ansatz, 

i^k{a,4>) = Ck{a)ipk{4>) ■ (36) 

We choose 

= e-'^^"' , (37) 

because real exponentials would lead to exponentially increasing wave functions for cf) —>■ ±oo 
that would not reflect classical behaviour. From (jHn|l one then gets the following equation 
for the Ck (primes denote derivatives with respect to a), 

a^C'^ + aC'k + ^{a^ -k^)Ck = . (38) 

Solutions of this equation are Bessel functions Zk/2h{0'^ f^h). However, we have to impose 

the boundary condition that ipia, (p) in order to reflect the behaviour of the classi- 
cal trajectories (fT^ which have in configuration space a minimum with respect to a. We 
therefore have to choose the Bessel function Jk/2h{o? /2h) with k > 0. 

The connection to the classical solution (fT3j) should be performed through a formal WKB 
limit '7i — > 0'. We thus have to look for an asymptotic expansion of J where both the 
argument and the index are large. We use the expression [s^l 
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and set v = k/2h, z = a'^/k. The choice for ( depends on whether > 1 oi < 1. Let us 
consider first the case z"^ >1 which corresponds to a^/A;^ > 1. From one sees that then 



2/3 



We also use the asymptotic expression for the Airy function occurring in (jHUj) . see js^, 

2/3 



Ai 

where 



2h 



C ~ TT 



-1/2 



2h) ^ 



sin^fc , (41) 



0. - -l-,e- + \ . (42) 

The classical trajectory is then recovered through the principle of constructive interference: 
We look for the extremum of the phase 

S,^e,±^ (43) 

of the total wave function with respect to k. One then easily finds that the requirement 
dSk/dk = at k = k leads to (fT3jl . One can thus identify Cp = k. 

What happens for z'^ < 17 As one can easily see from the corresponding expression in 
list - C < and the Airy function decays exponentially. This is as expected, since a'^/k'^ < 1 
corresponds to the classically forbidden region. 

One can also easily check that 5*^, Eq. PHj) . is a solution of the Hamilton-Jacobi equation 
arising from (jH) through the substitutions tt^ dSk/da and tt^ dSk/dcf). 

In the case of the conventional scalar field, one gets a change of sign for the /c^-term in 
fjHHjl . The solutions for Cfc(a) are then the Bessel functions Jik/2h{ci^ f^fi) and J_jfc/2n(a^/2/i). 
Since there are no classically forbidden regions, both solutions seem to be allowed. It can 
again easily be checked that the classical solution follows in the formal limit — > 0' 
from the principle of constructive interference: One gets the two branches of (fTH|) from the 
two Bessel functions. This suggests to use one or the other Bessel function if one wants to 
avoid interferences (and thus non-classical behaviour) at large a. Since (j30|) is hyperbolic 
for i = 1, one is free to impose boundary conditions at constant a, that is, one can either 
impose one packet or two packets there, depending on whether one wants one branch of the 
classical solution to be represented or both. 

In the phantom case discussed above, the Wheeler-DeWitt equation is elliptic; one there 
only imposes the boundary condition that ip goes to zero at a ^ and that it is at most 
oscillating at the other boundaries. This fixes the solution to be Jfc/2a(a^/2/?.) or a superposi- 
tion thereof. Explicitly, one would consider the following superposition for the construction 
of a wave packet, 

poo 

^P{a,^)= / dfc A(fc)e-*'*/V,/24aV2^) , (44) 
Jo 

where A{k) is a function of k that is peaked around a particular value k (e.g. a Gaussian). 
One would not expect the packet to exhibit dispersion near the minimum of the classical 



12 



trajectory, since the phase of the Bessel function is not rapidly varying with respect to k, 
in contrast to the case of a massive scalar field discussed in [35|. We shall, however, expect 
the occurrence of a dispersion at large values of a. We shall discuss this explicitly for the 
more realistic case in Sec. Ill C below. 

Making an analogy to ordinary quantum mechanics, one would compare the solution in 
the elliptic case to an 'initial wave function' iplt = 0, x), whereas the hyperbolic case would 
correspond to the time evolution tlj{t,x), since one would have for (j3U|) a distinguished set 
of foliations with respect to an intrinsic time defined by the scale factor. This intrinsic 
time could be used as a physical time with respect to which, for example, further degrees of 
freedom could be evolved, cf. Sec. IV. 



C. Quantum phantom cosmology — exponential phantom potential 

For non-zero, exponential potential as in Sec. Ill CI the Wheeler-DeWitt equation is 
most conveniently solved after a transformation to new variables in such a way that the 
potential cancels in front of \1'. This is obtained by first transforming to light-cone type 
coordinates zi ^ a + z\ = a — For £ = 1, these are just the characteristics of the 
Wheeler-DeWitt equation. The equation now takes the form 

from which a transformation to new variables can be made such that / (zi, Z2) is canceled. 
This corresponds to a transformation to the variables 



<P) = 7 — 72 cosh(X) + sinh(X) 



4>) = 7 — -r„ sinh(X) + i— cosh(X) 

3 \Vi ^/6 

where X = v^(30 — i^^a). For both the phantom and the ordinary field, U£ and ve are 
real. The Wheeler-DeWitt equation in these variables takes the simple form 

Making a WKB-approximation ansatz, \1' = Ce^^'^, one obtains at lowest order the 
Hamilton-Jacobi equation 



duf ) V dv. 



(47) 



This is solved via a separation ansatz by ^ofc = ku^ — \J l{k? — Of course, the Hamilton- 
Jacobi equation is also solved by actions carrying different signs in front of and vi. These 
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are obtained from the one chosen above by rotations in the (m^, f^)-plane. For i = —1, all 
solutions can be mapped onto each other in this way. This is an obvious consequence of the 
rotational symmetry of Eq. ()47|) for £ = —1. As Mi > (recall that A < for i = 1) for 
the conventional scalar field, here only two solutions can be mapped onto each other. 

From the classical action Sok, the equations of motion are obtained via ^1^1;^,=^ = c. 
(Note that 5*0^ evaluated ai k = k is always real.) For the special case c = and 



one obtains the classical trajectories 



6 



0(«) = i^a , (48) 

cf. (EHl). 

Plugging this lowest-order ansatz into the Wheeler-DeWitt equation, one finds that the 
equation is already satisfied exactly. Thus we get the following exact wave packet solution 
to the Wheeler-DeWitt equation, 

^{ue, ve) = j dkA{k) (^Cies(^'"^-V^(^''^) + C2e-^('=«*- V^C^'-i)-^)) . (49) 

By construction, the classical trajectories can be recovered from this equation through the 
principle of constructive interference. We choose for the amplitude a Gaussian with width 
a centered around A;, 

m 



Taking Ci = C2 for definiteness, one obtains wave packets of the form 



where a Taylor expansion of S^k has been carried out around k (primes denoting derivatives 
with respect to k) and the terms of the order {k — k)^ in the exponent have been neglected. 
(For simplicity, in this expression Sok{k) = Sq.) This can be done if the Gaussian is strongly 
peaked around k, that is, if a is sufficiently small. Since S'^i^ik) = gives the classical 
trajectory, the packet is peaked around it. For the conventional scalar field as well as for the 
phantom field, the wave packet thus follows the classical trajectory but spreads as vj — > 00. 
This can be recognized from ()50|). since the term proportional to [>S'q^(A;)]^ in the width of 
the Gaussian increases without limit, 

S'LCk) = , 3 • (51) 



{l{k^ - 1)) 



2 



It is even more obvious from the absolute square of the wave packet (neglecting for simplicity 
the complex conjugate part in (|^). 
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The spreading occurs due to the non-trivial dispersion relation, that is, due to the fact that 
5*0*; depends non-linearly on k. The semiclassical approximation is thus not valid throughout 
configuration space. 

For the phantom field we have U-i — > — oo, V-i oo when we approach the Big-Rip 
singularity. This singularity thus lies in a genuine quantum region. Since for £ = —1 one 
has 

it is obvious that the occurrence of the non-trivial potential is responsible for the dispersion. 

The Big-Rip singularity is thus 'smoothed out' — when the wave packets disperse, we 
can no longer use an approximate time parameter; time and the classical evolution come to 
an end, and one is just left with a stationary quantum state. This corresponds to quantum 
gravity effects at very large scales. Hitherto such a case has only be encountered near the 
turning point of a classically recollapsing universe, as a consequence of the demand that the 
wave function go to zero for large scale factor [35. 41|. 



Due to the fact that mi > for the conventional scalar field model, here two inequiv- 
alent actions exist. Apart from the wave packet constructed from the function Sok = 
kui — ~ 1^1, one gets a second wave packet constructed from Sok = —kui — i/fc^ — \vi. 
Moreover, the entire (a, 0)-plane is mapped into only one quarter of the f i)-plane. One 
would therefore require the wave packet to vanish at the boundary of the physical region. 
The only solution satisfying this requirement is naturally the trivial one. To get a non-trivial 
solution, one has to lessen the boundary condition and require = only at the origin of 
the fi)-plane. The fact that the wave packet does not vanish at the ui = and f i = 
line is due to the non-normalizability of the wave packet in both a and 0, which in turn has 
its origin in the fact that the classical trajectory has no turning point. 

The implementation of this condition results in a wave packet which vanishes at the Big- 
Bang singularity, \E' — > as a — > —oo, and spreads for large a. The Big-Bang singularity does 
therefore not exist in the quantum theory. In the phantom field model, no such restriction 
occurs due to the fact that the entire f_i)-plane represents the entire (a, 0)-plane. The 
wave packets for both the phantom and the ordinary scalar field are depicted in Figure 3. 



D. Quantum phantom cosmology — scalar field fluid and negative cosmological 
constant 

For the model discussed in Sec. Ill D[ the classical solutions require a potential of the 
form V{(j)) = Vq cosh^(0/F), cf. (j^ . The Wheeler-DeWitt equation therefore reads 

The classical singularities lie in a region of large In order to study the quantum behaviour 
there, it is thus sufficient to approximate the potential for large |0|, 

V-(^) « ^e^ , (54) 
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V 



FIG. 3: The amplitude of the wave packet ()5U() for an exponential potential solution of the WDW 
equation for ^ = 1 (left) and I = —1 (right). Here, h was set to unity and parameters a = 0.1 and 
X = k/2 have been chosen. The wave packet for the phantom field model is seen to spread near 
the classical singularity. For the scalar field model one has ^' ^ at the origin. In each sector 
corresponding to one copy of the (a, cj)) plane, the same wave packet propagates. 



where in the following the upper sign refers to positive 0, and the lower sign to negative 
0. This makes the problem very similar to the one of Sec. IIII CI The Wheeler-DeWitt 
equation is here simplified by a transformation on the variables 



U£{a,(f)) = — -= ^ I cosh (X) =1= ^ sinh (X) 



3^2 1 - ^\ 3FVi 



veiaA) - (^-^ sinh (X)±- cosh (X) 

where X = y/i (30 ± ia/F). In these variables, we recover the form 

Again, one obtains a solution from a WKB ansatz. The Hamilton- Jacobi equation is again 
given by ()47p (this equivalence is, of course, only formal, since ui and ve are defined differ- 
ently). It is again solved by Sok = ku^ — iik"^ — l)f^, where the remarks of Sec. IIII CI 
concerning the choice of action apply here as well. The equations of motion obtained for 

Mofel - 
dk \k=k 



are 



£ D 

0(a) = V ^ ^^'^ ■ ^^^^ 

This solution coincides approximately with the classical solutions (|28p: If one approximates 
for £ = —1 for large a, one gets (± label the different branches of the classical solution) 
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= ± 




a±F\n{2A) , 



(57) 



where a > 0. Therefore, the hmit of large positive is obtained on the 0_|_-branch, while the 
limit for large negative is reached on the 0_-branch. On the other hand, an appproximation 
for small a in the case i = 1 yields 



where a < 0. Due to this, the limit of large positive is obtained on the 0_-branch, and for 
large negative on the 0_|_-branch. Thus the solution to the approximated Hamilton-Jacobi 
equation (jlTjl coincides with the approximation of equation (j2HI)- Of course, a special choice 
for k has to be made to fix the onset. The fact that for £ = — 1 large correspond to large 
a, and for i = 1 large correspond to small a is due to phantom-scalar field duality. 

With the help of the classical action 5*0^, the approximate Wheeler-DeWitt equation can 
be solved. Again, the WKB ansatz satisfies the equation exactly. The wave packet is of the 
same form as in Sec. IIII C\ with a different definition of ue and ve and another choice of 
the center of the Gaussian, k. As in the case of vanishing cosmological constant, the wave 
packet spreads for f| — > cxd. The Big- Rip singularity in these variables occurs at v'^^ — »• oo, 
U-i —>■ oo. Thus, again, the singularity is hidden in a quantum regime and the semiclassical 
approximation is not valid throughout configuration space. 

Due to the restriction D < 3 for the i = 1 model, the same remarks concerning the range of 
the coordinates as in Sec. EEl apply here. So at the Big Bang. * 0. I„ analogy to Q 
one would expect quantum effects to occur also in the region of the classical turning point. 
This will be addressed in a future publication. 

IV. DISCUSSION AND OUTLOOK 

In our paper we have applied the formalism of standard quantum cosmology (using the 
Wheeler-DeWitt equation) to a situation where phantom fields are present. This is of 
interest because there are novel features with regard to both the structure of the equation 
(elliptic or ultrahyperbolic instead of hyperbolic) as well as the presence of new scenarios 
(Big-Rip singularity at large scale factors in the classical model). In fact, one of the most 
intriguing features is the possible occurrence of quantum effects for large scale factors. 

For various models we have determined and discussed the classical trajectories in con- 
figuration space. We have then considered the corresponding Wheeler-DeWitt equations; 
we have given various solutions and addressed the classical limit as well as the behaviour of 
wave packets following the classical trajectories in configuration space. We have found that 
the packets disperse in the region of the classical Big-Rip singularity. This singularity is 
thus 'smeared out' by quantum effects at large scale factor. Once the wave packets disperse, 
no approximate time parameter can be defined [s^ and the classical evolution terminates 
in a singularity-free way. 




(58) 
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For the conventional scalar field model we have found that the wave packet vanishes at the 
Big-Bang singularity due to the implementation of appropriate boundary conditions. In this 
way, the Big-Bang singularity is removed from the quantum theory. This is similar to the 
avoidance of the singularity in models of loop quantum cosmology [3l| and shell collapse js^ . 
Without this boundary condition the wave packet would just have approached the region 
a — oo without spreading; this lack of dispersion is a result of the Wheeler-DeWitt 
equation taking the form of a free wave equation in this limit. 

The present work can be extended in various directions. The next step would be to add 
other 'conventional' scalar fields and to investigate the full quantur n dy namics. In particular, 
this would be of importance for a discussion of the arrow of time |23]. In order to define an 
appropriate entropy it is necessary to introduce a set of inhomogeneous degrees of freedom. 
In the case of a classically recollapsing universe it has been found that the arrow of time 
is correlated with the scale factor of the universe, that is, the arrow of time must formally 



reverse at the maximal expansion j3^|4l|. (The reversal is formal because quantum effects 



near the classical turning point do not allow classical observers to survive this region.) It is 
of interest to investigate whether a similar behaviour occurs here. Regarding the classical 
evolution of the phantom fields depicted on the left-hand sides of the Figures 1 and 2, one 
would again expect a correlation of the entropy with increasing scale factor, that is, there 
would be no collapse followed by an expansion but only two separate branches of expansion 
separated by a quantum region. The classical evolution would then start out of the quantum 
phase at small scales and describe an expansion of the universe ended by a quantum phase 
near the 'Big- Rip region'. We shall address this scenario in a forthcoming paper. This will 
also deal with a possible quantum phase near weak singularities (sudden future singularity, 
generalized sudden future singularity, type III and type IV) as mentioned in the Introduction. 

We have based our discussion on the Wheeler-DeWitt equation of quantum geometro- 
dynamics. More recently, an alternative formulation of canonical quantum gravity called 
loop quantum gravity has gained considerable attention, cf. js^. A major prediction of this 
approach is the presence of a discrete structure for geometric operators. This formalism 
was applied to cosmology where it led to new features j3l|. Instead of the usual Wheeler- 
DeWitt equation one gets a difference equation for the scale factor. While near the Planck 
scale this equation gives different results from the (differential) Wheeler-DeWitt equation 
(and therefore can prevent the occurrence of the classical singularities), it coincides with it 
for higher values of a. It thus seems that near the classical Big-Rip singularity the same 
scenario emerges that has been discussed in the present paper. However, it is of interest 
to investigate quantitatively the differences and similarities of ordinary quantum phantom 
cosmology and loop quantum phantom cosmology. A first paper in this direction has studied 
the effective dynamics from loop quantum cosmology and its consequences . The results 
indicate that the Big Rip can be avoided. We hope to return to these and related issues in 
a future publication. 
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